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Abstract. This paper presents a methodology to stabilize some kind of Non- 
linear Control system known as Driftless, utilizing the concept of Pseudo- 
Kinetic Energy introduced in this work. Once this controller is applied to the 
Unicycle-type robot, stability is guaranteed with the salient property that the 
structure of the controller allows to solve in closed-form the trajectories of the 
vehicle. 

While the proposed controller only ensures stability (not asymptotic sta- 
bility) the obtained closed-form solutions will show a path to obtain in closed 
form the solutions for the general control problem of the unicycle. 

Some conclusions and future directions for research are also depicted. 



1. Introduction 

The stabilization of Driftless affine systems is a very important topic both in 
practical cases and theoretical analysis (Non holonomic vehicles, robotic systems 
in chain form,etc), because of that, they have received a lot of attention in past 
years (see for instance [5] and [2] with the references trough them), 

While all the available techniques for stabilization of Nonlinear Driftless systems 
consider time varying open-loop control laws ([3], [2]) or smooth controllers using 
Lyapunov Functions (see [?]); up to the authors knowledge the implementation of 
a controller for Driftless systems with the possibility of Closed-Form solutions of 
stable trajectories was never done in the past. 

In this way, the goal in this paper will be to derive a general continuous smooth 
stabilizing control law for Driftless systems, introducing the concept of Pseudo- 
Kinetic Energy. This control law will be applied to the well known Unicycle-like 
robot, yielding Closed-Form solutions for the trajectories of the wheel baseline 
center coordinates and the attitude angle via the solution of a second order Bessel's 
equation. 

Once the closed-form solutions are obtained, the impossibility of an asymptotic 
behavior for one of the wheel baseline center coordinates will be established in 
accordance with the Brockett's theorem (see [IT]). However, as it will be shown, if 
the same controller is used in such a way that the attitude is unstable (the robot 
continues rotating), then the wheel baseline center will be in fact asymptotically 
stable. 

This conclusion will guide an extension for this controller via the addition of a 
dissipative term which produces some sort of "shaking" behavior of the attitude, 
endowing the entire Kinematic Model of the vehicle with asymptotic stability. It 
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turns out that the extended controller is time-varying preventing collision with 
Brockett's condition. 

It is worth mentioning that controllability is not a condition to be checked in the 
application of this methodology, even though the controllable nature of the system 
is not discussed in this work, and the contrast with cases where the Lie Algebra 
based test fails while the system remains controllable is not focused in this paper. 

This paper is organized as follows: Section [2] introduces formally the class of sys- 
tems considered in the paper, Section presents general theorems for stabilization 
of Driftless systems, Section 2] apply the previous theory to the Unicycle and finally 
Section ?? extends the analysis for asymptotic stability of the Unicycle. 

2. The Problem Considered and Motivating Ideas 

In this section, the problem considered in this paper and some preliminary ideas 
related with the stability criterion used to design stabilizing smooth controllers will 
be established. 

The kind of Kinematic Models considered in this paper are the well known 
"Driftless" systems: 

(2.0.1) q(t) = S(q) ■ u(q), 

where q{t) £ 3?™ xl are the state space variables, u(q) S 3? fcxl is the admissible 
control set which is at least C 1 (M), where M could be any manifold configuration 
for the state space q(t). 

The challenge is to develop a systematic technique with two main characteristics: 

• Stabilization of the system (|2.0.1j) in Lyapunov sense. 

• The possibility of addition of dissipative terms in order to make the system 
asymptotically stable. 

The first item will be reached introducing the concept of Pseudo-Kinetic Energy 
and showing that such a definition allows obtaining a Lyapunov function which 
ensures stability for Driftless systems. Realize at this point that asymptotic stability 
is not yet addressed. 

A special Lyapunov function is obtained by introducing the concept of Pseudo- 
Kinetic Energy for a subclass of Driftless systems. It is showed that the Unicycle 
is in this category. Moreover, the technique applied to the Unicycle yields a stable 
system and also allows to have closed-form solutions for the trajectories of such a 
system. 

Analyzing these closed-form solutions is concluded the impossibility of asymp- 
totic stability with this controller. The problem could be solved with the inclusion 
of a dissipative term. In fact, by defining a general dissipative term for the attitude 
angle of the Unicycle, it will be shown in an analytic way that the Kinematic Model 
will turn to be asymptotically stable. 

3. A Stabilizing Controller for Driftless Systems 

From now on, stability will refers to stability in the sense of Lyapunov. So, a 
stabilizing controller should satisfy the Lyapunov theorem which yields a sufficient 
condition for a Nonlinear system to be stable (see [10], pp. 188-189). Realize that 
the concept of 12-norm plays a fundamental role because it is a scalar quantity that 
allows to measure how near is the system to a stable equilibrium point. 
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In this way, since the focus is to drive the system to an equilibrium point (q = 0), 
and recalling that the Barbalat's lemma predicts that if an integrand is going to 
zero anytime the integral is bounded, and with afore mentioned the idea of the 
12-norm in mind; it is reasonable to consider the following integral 

(3.0.2) / m)M)) ■<% = Bounded, 

Jo 

where (., .) means the inner product and (<?(£)> <K£)) = \\q(t)\\ 2 . 

The machinery in what follows will derive a methodology to design a Control 
law u(q) for the Driftless Kinematic Model introduced in equation (|2.0.ip . In this 
way, the first result is an extension of the well known Barbalat's lemma - all the 
proofs are given in the appendix: 

Lemma 3.0.1. Let if : 5ft — > 5R be a positive definite junction. Suppose that 
t 

linit^oo J <p(t) ■ dr exits and is finite. Then, ip(t) — > as t — > oo. 
o 

This lemma shows that one way to guarantee q(t) — ► when t — ► +oo is to 
ensure that the integral in (|3.0.2p when t — > +oo is bounded and this motivates the 
following result: 

Lemma 3.0.2 (Stability Criterion). Given a Control System of the form q(t) = 
S(q(t)) ■ u(q), where S(.,.) is a smooth collection of vector fields which belong to 
sjjnxfe yjnfo q — o an equilibrium point. Then, this system is stable if and only if the 
quantity E* = J Q (<?(£),<?(£)) ■ d£ called Pseudo-Kinetic Energy is bounded. 

Next section is using this Lemma for a subclass of the more general Driftless 
one, providing a systematic methodology to obtain stabilizing controllers for such 
a systems. 

4. Stabilizing the Unicycle 

Based on the Lemma 13.0.21 in previous section, the main result in this paper for 
a special subclass of driftless systems is presented: 

Theorem 4.0.3 (Controller Calculations). Given a Driftless Control system to the 
form R2.0.1\) , such that \\Si\\ =1, % = 1,2, . . . ,fc and Si ■ Sj = 0, i ^ j where \\.\\ 
is the l2-norm - see [15 for details on matrix norms- and S = [Si, S%, . . . , Ski- 
Then the controller Ui{q) = p ■ q(t)' - Si, p < 0, stabilize the system q(t) = 
S{q{t))-u{q). 

Applying this result to the Unicycle, it will be showed that this kind of controller 
allows to find Closed-Form solutions for the trajectories of the system where a 
posterior modification via " dissipative" terms will produce asymptotic stability. 

4.1. Closed Form Trajectories. Considering the Kinematic Model of the Unicycle- 
see [10] pp.529: 



Also consider the ideas in the paper [8] 
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(4.1.1) 



cos(0) 
sin(0) 
1 



S(q) 



where q(t) = [ x c (t) y c (t) 8(t) 1 with X = [ x c (t) y c (t) 1 the position of 
the wheel baseline center of the robot and 0(t) the attitude of it. 
Clearly this system has the property: 



'Si = [cos(0) sin(0) 0]' 
S 2 = [0 1]' 

l|Si|| = l 

11^11=1 



Then applying the controller in Theorem 14.0.31 we have: 



cos(6») 
sin(6») 



where p < 0: 



(4.1.2) 



u 


(«) = 


p.X>- 






P 


x c (t) 




cos(0) 


Vc{t) 


= P ■ 


sin(6>) 



(cos(0) ■ xc + sm(6) ■ yc) 



[6(t) = p-6(t) 



Realize that the dynamic of the attitude (0(t)) can be decoupled from the rest 
and it is clearly asymptotically stable (globally) to zero. In this way, the con- 
cerned in what follows will be the position of the wheel baseline center of the 
robot (x c {t), y c (t)): 



x c (t) 

Ut). 

X(t) 



cos(0) 
sin(0) 



(cos(0) • xc + sin(0) • yc) 



This subsystem can be transformed conveniently via a state space transformation 
as follows: 



(4.1.3) 



x c (t) 




cos(0)" 


Vcit). 


= p ■ 


sin(0) 



[cos(0) sin( 



x c (t) 

ydt) 



X(t) Y(9) 

Moreover, it is convenient to write Y(0) as follows: 



Y(9) = 



cos((9) - sm(6) 
sin(6») cos(0) 

R(6) 




■R(ey. 
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It is important to notice that R(9(t)) is a rotation matrix about 9 and because 
of that is always invertible and smooth. In this way let's consider the following 
change of coordinates: 

Definition 4.1.1 (Change of Coordinates). X(t) = R{6) ■ Z(t). 

Notice that this transformation is a Lyapunov transformation (see [11] . pp. 107- 
108 for details) and because of that preserves the stability properties. 
Applying this transformation to (|4.1.3| with the definition above: 

(4.1.4) x(t) = P ■ R(ey • a • R(ey ■ x(t) 

Z{t) 

On the other hand, it is possible to derive an expression for Z{t): 

Z(t) = R{9)' ■ X(t) + R(e)' ■ X(t) 
Replacing X(t) from equation (|4. 1 .4j) : 



Notice that: 



Z{t) = R{6)' ■ X(t) + p- A ■ Z{q) 

R(9)' = -e\t)-\-^ cos( ^ 

w w — cos(fl) — sm(t>) 



R(8- 



then: 



Z(t) = -6(t) ■ R{0 - -)' ■ R(9) ■ Z(t) +p-A- Z{q) <^ 

z(t) = -[6(t)-R(e-~y-R(e)- P -A]-z(t) 

Using the multiplicative property on rotation matrices: 



this leads: 



Also: 



R{e-\y = RC-)-R(-6). 



Z(t) = -[9(t) ■ R(-) ■ R(-9) ■ R{9) - p ■ A] ■ Z(t) 



R(-6) ■ R{9) 



1 
1 



then: 



Z(t) = -[9(t)-R(-)-p-A]-Z(t). 



This can be recast follows: 
Z(t) = -[9(t) 



' -1 " 




' 1 


" 


1 


- p ■ 









Z(t) 
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In summary: 



(4.1.5) 



X(t) = p ■ R(9) ■ A • R(0)' ■ X(t) 

-p -m 

9(t) 
zi{t) z 2 {t) 



Z(t) = 
Z(t) = 



Z(t) 



This first order set of Ordinary Differential Equations (ODE) can be decoupled 
into two second order ODE's as follows: 




9{tf)- Zl {t) 



e(t)\ ■ (f \ , , n e(t) 

z 1 (t)-p-z 1 (t) 

<?(*) 

Moreover, with the Kinematic Model for 9{t) in equation (|4.1.2p , this set of 
ODE's yields: 

(z 1 (t) = 2- P - h(t) - P 2 -{i + eit) 2 ) ■ z x {t) 

(4.1.6) | z 2 (t) = Zllt) e ^ Zlit) 

[e(t) = p- e(t), P <o 

In the Appendix it is shown that (|4. 1 .6|) can be converted to a second order 
Bessel's equation which posses general solutions J p (x) and Y p (x) called Bessel's 
functions of first and second kind respectively -see [12] . pp. 228-230. 

From now on and for the sake of clarity, it will used p = — 1 - otherwise the time 
change: r = p ■ t transform the problem to this form. In this way and using the 
solutions in the Appendix, the closed- form expressions for z\ (t) and z 2 (t) using the 
Kinematic Model for 9 in equation (|4.1.2|) . arise: 



'z 1 (t)=8-[C 1 -Jo(8) + C 2 -Y (9)} 



(4.1.7) ^ o = P • e 

_ ii(t)+zi(t) 
6(t) 

where C\ , C 2 are arbitrary constants depending on the initial conditions. 



z 2 (t) 



4.2. Asymptotic Analysis and Further Conclusions. The analysis of the be- 
havior of the closed-form solutions obtained when t goes to infinity will clarify 
some issues regarding asymptotic stability for the different immersed submanifolds 
forming the orbits. 

In this way and using the formulas for small arguments presented in [12] . pp. 
248-249, the analysis for zi(t) from equation (I4.1.7|) : 



(4.2.1) 



lim zi (t) = lim 8 ■ [d ■ J (9) + C 2 ■ Y Q (8)} = 



8^0 



providing that Jq(9) is bounded V£? € 3? and lim^o 9 ■ Yq(9) = 0. 

The first conclusion is that z\(t) is asymptotically stable, however as it will be 
shown, this is not the case for z 2 (t) . In this way, and using the property for Bessel's 
functions (see [12] . pp. 199 and pp. 226): 
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\lE J n(x) = 2 t-Jn{x)- J(n+l){x) 
\&M X ) = 3 ■ ( Y (n-l)(x) ~ Y [n+l) {x)) 

In the case of equation (|4.1.7p , a — 1 and n = 0, 2!a(t) yields: 



(4.2.2) 



a*(t) 



C 2 



■ [-ci ■ J ( i) (0) + -y • (y ( -D(0) - 

The formulas for small arguments -see [12], pp. 248-249, lead: 



JlimeWnW = (f)" • r^ry 
\u m ^ o Y„(0) =-(§)«. EM 

where T (.) is the Gamma function. Then if 9 — > 0, in equation (|4.2.2p : 



lim zofi) = lim \—Ci ■ r-r — Co • 

e^o v 7 e^o 1 2-T(2) 



6> 2 -r(-i) 

4 • 7T 



+ C 2 - 



r(i) 



Since T(l) = T(2) = 1 and is finite, then: 



lim z 2 (t) = C 2 ■ - 

To investigate the possibility for Ci = endowing Z2 with asymptotic stability, 
Definition 14. 1 . 1 1 is used, taking into account that C\ and Ci are arbitrary constants 
depending on the initial conditions as pointed out in equation (|4.1.7p : 

X(t = 0) = R(9(0))- 

' 9(0) ■ J O (0(O)) 9(0)-Y Q (9(0)) 
' -9(0) ■ J (9(0)) 9(0) ■ i ■ (Y^(9(0)) ~ n(«(0)))J 

where the closed-form Kinematic Model for z 2 form equation (|4.2.2[) was used. 
In this way, if C2 — is settled: 

X(t = 0) = 9(0) ■ Ci • J ■ R(0(0)) 

which is not allowing any initial condition as required to have 3? 2 a the set of 
attractors -see [13] for a complete discussion on attractors. This fact agree with the 
Brockett's condition regarding the impossibility to stabilize the unicycle by means 
of continuous controllers -see [T7] . 

A last important conclusion arises if the set of attractors is considered: 



X(t->oo) =R(0)- 
1 

That means that the set of attractors agree with the whole y c axis. On the other 
hand, if p = 1, then applying the formulas for 9 — > 00 , it can be seen that z\ and 
Z2 tend asymptotically to zero leading instability in 9(t). 

Next section utilizes this conclusion, extending this controller adding a dissipa- 
tive term, endowing the system with an asymptotic behavior. 
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5. Comparison with other Methodologies 



Once the methodologies in this paper have been presented, it is important to 
compare them with some of the available in the literature. In fact, the survey in 
PQ will be followed in order to classify them. 

5.1. Brockett's Theorem. As pointed out by Brockett in his theorem in [IT) , 
systems like the Unicycle can not be asymptotic stabilized using smooth controllers, 
however even when it is simple to arrive to this conclusion using the Brockett's 
theorem it is not simple to expose in a clear way the impact of this condition for 
every particular controller designed for the Unicycle. 

In this sense the methodologies in this paper shows clearly the impossibility in 
using smooth controllers in order to produce asymptotic stabilization for the unicy- 
cle since closed-form solutions are available and moreover the asymptotic analysis 
of previous section determine this condition in an alternative way. 

5.2. Comparison with Smooth time-varying techniques. While Method 1 
and 2 in [I], pp. 5-13 need the evaluation of Lyapunov functions to construct sta- 
bilizing controllers for Driftless systems, the techniques in this paper is systematic 
and is no based in the Lyapunov theory. However when compared with the method- 
ology of Pomet-Coron -see [I], pp. 14- the method of the present paper serves as a 
start point for the three step algorithm in pQ, pp. 16, specially step 1. 

5.3. Comparison with Homogeneous time-varying Feedback. Even when 
the theory under this subject is well developed and understood -see p], pp. 20-28- 
the theorems in the present paper confirm some results developed using Homoge- 
neous time-varying Feedback, in particular compare Theorem 10 in [T], pp. 20 and 
Condition ?? in the present paper. 

5.4. Comparison with the Transverse functions technique. This technique 
is similar to the heart of the methodologies presented in this paper, in fact the 
stability Theorem 14 . . 31 can be thought as a Transverse function or in other words, 
can be depicted as the control system running perpendicular to some Submanifolds 
just like the Transverse function technique. 

In this sense the Transverse function technique is more flexible and better devel- 
oped but the present paper addresses a closed-form expression allowing a further 
analysis of the stability properties for the Unicycle. 

Clearly the lack of flexibility omnipresent in the technique of the present paper 
is a disadvantage. 



A slightly different model than the one consider so far can be also integrated to 
obtain closed- form solutions, let's consider the case presented in [5] and write the 
control vector as: u\,a ■ u%, where the parameter a is the distance from (x,y) to 
the center wheel base-line. 

The, utilizing the same procedure as before: 



6. The complete solution 



X(t) = R(9) ■ [I - R(-) - R(6(t)) ■ R(-6(0)) ■ R(-)} ■ R(0(0)Y ' + 
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where X(t) = [x(t),y(t)]' , u(t) — [ui,a ■ U2]' 
analyze the effect of the parameter a — > if limt- 



. In this way, it is possible to 
,0(t) = : 



Umt^ooX(t) 



1 -1 
1 1 



POO 




'l -1" 


Jo 




1 1 



R(-d(0)) ■ R(- t 



R(-e(a)) ■ R(-) 



■X(0) + 
u{a) ■ da 



Clearly, when a — one control input is lost in rtdert to steer the vectroX to 
zero. This explains why the Brockett's condition is working. It is also worth to 
notice that this case is of more general application (since a = include the previous 
case non-holonomic robots) and any control strategy like model predictive control, 
for instance , can be applied. 

7. Conclusions and Future Work 

A methodology to stabilize Driftless Nonlinear Control systems was presented, 
by applying this controller to the well known Unicycle-like vehicles, closed-form 
solutions for the trajectories were obtained via the solution of the second order 
Bessel's ODE. 

The analysis of the closed- form solutions showed that the coordinate x c is asymp- 
totically stable while y c does not. However, if the simple stabilizing law for 9: 
9{t) = p ■ 6(t), p < 0, which arises with this technique is considered in the op- 
posite side, i.e. p > 0, then both coordinates x c , y c turn to be asymptotically 
stable. 

This conclusion can be exploited in tow main ways: 

• Steer the vehicle using p > until the coordinates x c , y c are closed enough 
to the origin and then to switch to p < stabilizing the vehicle. 

• Adding a dissipative term making the controller time- varying such that the 
closed-loop system is asymptotically stable 

As future work, several directions can be depicted, in fact, since the same con- 
stant p is used in both components of the input control vector for the Unicycle, it 
would be of interest to relax this study utilizing two different constants. On the 
other hand, this closed-form solutions could be exploited for formations of robots 
where the incorporation of constraints should be considered. 



8. Appendix 

This appendix contains the proof of all the results along the paper. 
Proof. (Lemma \3.0.1\) 

t 

The first step is to write lim^oo J f{r) ■ dr as an infinite sum: 



t „ i-T 



\im J <f{r) ■ dr — J <f{r) ■ dr, 



—^(i-l)-T 

where T ^ 0. As it is well known for infinite series -see [IB] , a necessary condition 
for convergence is the general term tending to zero: 
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■T 



lim / <p{t) ■ dr = 0, 

^oo J 



l — >oc 

(i-l)-T 

Since T is finite and ip is positive, then the only possibility is (p(t) — > as t — > oo 
(or i — > oo). This completes the proof. □ 

Proof. (Lemma V3.0.2\) 

Let first consider that the quantity E*(oo) — f£°{q(£), ?(£)) ' i s bounded 
and prove that the system is stable. To do this, it is required to prove that the 
function E*(t) = J (<?(£), <K£)) ■ d£ is strictly increasing and this can be done with 
the following useful result: 

Given a function f : Sft — > 5R + , then, cp(t) = f f(r) - dr is strictly increas- 
ing. 

The proof of this assertion is carried out by showing that f* 1 /(r) -dr < J o * 2 /(r) • 

dr if ti < t%. In fact, taking into account f* 2 f{f) ■ dr > 0, then J"' 2 /(r) • dr + 

/ tl /(r).dr>/ * 1 /(r).dr. 

Then, the following Lyapunov function is chosen: 

- 1 1 

v[q ' t} -EW) m^-y 

Realize that this function classify as a Lyapunov candidate function -see |10j . 
pp. 188-189. In this way, we have: 

' v(q = 0, t) = w(0, oo) = 0, Since q — is an equilibrium point. 
v(q,t) is decrescent 

v(q,t) >o 

v{q,t)=-^ w -E* c {t). 

Since E*(t) > then —v(q,t) > which completes the proof for this part. The 
counterpart of the theorem is obvious taking into consideration that only smooth 
vector fields are considered. □ 

Proof. (Theorem \l0l\j 

The proof is based in Lemma [3.0.21 showing that the Pseudo-Kinetic Energy, E*(t) 
is bounded. 

Previous the calculation of E*(t), it is important to notice: 
Ui{q) = -p - q 1 ■ Si 

i = 1,2, . . . ,fc k 
S' z -S,=0 : i?j ^q(t)=p-(J2S' i -S i )-q(t) 



S i 'Si — 1 
Then, E*(t) arise: 
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E* c {t) = p- I g(0 , -(E5 i .u i ).de = p.(||?(t)|| a -||g(0)|| 2 ) 

which is bounded provided p < 0. □ 

8.1. Transforming the second order ODE in Section l4.1l equation (|4.1.6[) . 

The focus in this digression is to transform the second order ODE in equation 
(|4.ip into a Bessel's ODE. To do this, the following transformation inspired in the 
technique presented in [12], pp. 229-230, is considered: 

Zl (t) = x(t) a ■ l(x(t)) 

Then, z\ and Z\(t) lead: 

'ix(t) = [a ■ xit) 01 - 1 ■ l(x(t)) + x(t) a ■ l'(x(t))} ■ x(t) 
z\{t) = [a ■ xit) *- 1 ■ l(x(t)) + x{t) a ■ l'{x{t))] ■ x(t)+ 
* +[a -(a -I)- x{t) a - 2 ■ l(x(t)) + 2 • a ■ x^- 1 ■ l'(x(t))+ 

^+x(tr-i"(x(t))}-x(t? 

where = l'(x(t)) and — l"( x (t))- Replacing this new variables into 

the second order ODE at (|4.1.6[) : 



x(t) ■ [a ■ xit) 01 - 1 ■ l(x(t)) + x(t) a ■ l'(x(t))}+ 

+x(t) 2 ■ [a ■ (a - 1) ■ x{t) a - 2 ■ l(x(t))+ 

2 • a ■ x^- 1 ■ l'(x{t)) + x(t) a ■ l"(x(t))} = 

2 • p • x(t) ■ [a ■ xit)*- 1 ■ l(x(t)) + x(t) a ■ l'(x(t)]+ 

~p 2 -(l + 9(t) 2 )-x a -l{x(t)) 

Reordering: 



x(t) 2 ■ x(t) a ■ l"{x(t)) + l'(x(t)) ■ [x(t) ■ x{t) a + 
a ■ x(t) 2 ■ xit)"- 1 - 2 • p • x(t) ■ x(t) a ] + 
+l(x(t))-[x(t)-a-x(t) a - 1 + 

+x(t) 2 ■ a ■ (a - 1) • x{t) a - 2 - 2 • p • x(t) ■ a ■ x(t) a - 1 + 
+p 2 ■{l + 9{t) 2 )-x(t) a } =0 

In order to transform this equation into a second order Bessel's ODE to the 
form x 2 ■ d Jffi + x ' ^TtT + ( a;2 — n2 ) 'K x ) = 0, the following equalities should hold: 
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.1.1) 



.1.2) 



x(t) 2 ■ x(t) a 



x{t) 2 



.1.3) 



/(*(*)) 

x[t) ■ x(t) a + 2 • a ■ x(t) 2 ■ x(t) a - 1 + 

—2-p - x(t) ■ x(t) a — T T^TTyT 

x(t) ■ a ■ xit)"- 1 + x(t) 2 ■ a ■ (a - 1) • x(t)"- 2 + 
-2- p-x(t) ■a-x(t) a - 1 + 

x(t) 2 - n 2 



V ■ii + e{ty)-x{t) a = 



where f(x(t)) is any smooth no-null real function Vt £ 3?. In order to solve this 
algebraic system for a, n, f(x(t)) and x(t), equations (|8.1.2p and (|8.1.3[) are used 
equaling for the term x(t) ■ x{t) a : 



2 ' a ' /(i(t)) 2 ' P ' ^f*) ' X W) C /'!.- ( / I > 

= 5 9ggr 2 - 2 ■ p ■*(*)■ *(*)"+ 



V • (1 + 0(£) 2 ) • a;(i) Q • ^ - (x(t) ~" } • ^ 



f(x(t)) ' a 



that is: 



1 " <<^ Q ^ ■ /(*(*)) - (a;(t) 2 - n 2 ) • ^ 



a-x(t) = p 2 ■ (1 + 0(t) 2 ) ■ x(t) 



Guessing a possible solution for x(t): x(t) = 9{t) = 9 ■ e p t with p = ±1: 
fa = A- /(*(*)) +n 2 -i 

Replacing that function f{x(t)) into equation (|8.1.1[) : 



(8.1.4) 

Using equation (|8.1.2[) 



±(t) 2 ■ x(t) a — f( x (t)) 

x{t) = 6{t) = 9 • e ±4 



=>■ a — n = 1 



.1.5) 



1 + 2- a 



( a _ si) . a 



Replacing (j8T4|) into ([8713]) : 



a = 1 

n = ±V« 2 -1 = 



In summary: 
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( p = ±l 

i x(t) = 9(t) = 6»o • eP 4 
; a = 1 
' n = 

Zi(£) = x(t) a -l(x(t)) 

x ..i^ +x .*M> +[x 2_ n . ) . l{x) = Q 
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